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We study the formation of bright solitons in the impurity component of Bose-Einstein condensate- 
impurity mixture by using the time-dependent Hartree-Fock-Bogoliubov theory. While we assume 
the boson-boson and impurity-boson interactions to be effectively repulsive, their character can be 
changed spontaneously from repulsive to attractive in the presence of strong anomalous correlations. 

In such a regime the impurity component becomes a system of effectively attractive atoms leading 
automatically to the generation of bright solitons. We find that this soliton decays at higher tem¬ 
peratures due to the dissipation induced by the impurity-host and host-host interactions. We show 
that after a sudden increase of the impurity-boson strength a train of bright solitons is produced 
and this can be interpreted in terms of the modulational instability of the time-dependent impurity 
wave function. 

PACS numbers: 05.30.Jp, 67.85.Hj, 67.85.Be 


I. INTRODUCTION 

A soliton is a self-focusing solitary wave that main¬ 
tains its shape while it travels at constant speed and 
arises from a balance between nonlinear and dispersive 
effects. Bose-Einstein condensates (BECs) constitute a 
best environment for studying nonlinear macroscopic ex¬ 
citations in quantum systems. Excitations in the form of 
dark solitons with repulsive interaction have been real¬ 
ized in [in one and half decade ago. 

Bright solitons have been observed in BECs of ^Li in 
quasi-one-dimensional (ID) regime Si- The obser¬ 
vation of bright solitons was therefore possible only by 
means of the Feshbach resonance and then tuning of the 
interactions from repulsive to attractive during the ex¬ 
periments. In the experiment of Strecker eto/.i, the for¬ 
mation of the bright soliton trains has been interpreted 
as due to quantum mechanical phase fluctuations of the 
bosonic field operator Q. Bright soliton trains can be 
also generated in a BEC embedded in a quantum de¬ 
generate Fermi gasQ as a result of a competition be¬ 
tween two interparticle interactions: boson-boson colli¬ 
sions which are effectively repulsive and boson-fermion 
collisions which are attractive. 

In this paper we propose a novel scheme to realize 
bright solitons in quasi-ID atomic quantum gases. In 
particular, we study the formation of bright solitons in 
the impurity component of BEC-impurity mixture at fi¬ 
nite temperature by employing a versatile model known 
as time dependent Hartree-Fock-Bogoliubov (TDHFB) 
[ 2 ^ [ 23 . Experimentally, such mixtures have been al¬ 
ready realized with a medium composed of either bosonic 
i-El or fermionic atoms [iM|. Impurities in a Bose 
gas (Bose polarons) have been the subject of intense the¬ 
oretical studies. An important feature of these 
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mixtures is that when neutral impurity atoms immersed 
in a BEC can spontaneously form a self-localize state. 
This localized state, within the strong coupling approach, 
exhibits a solitonic behavior at both zero and finite tem¬ 
peratures [ 13 , [15 in quasi-ID geometry. These solitons 
are reminiscent of the well known optical wave solitons 


Here, we have pointed out that the bright solitons 
can be created spontaneously in the impurity compo¬ 
nent for high anomalous density. This latter quantifies 
the correlations between pairs of condensed atoms with 
pairs of non condensed atoms and rises with interactions 
strength. It was shown that the anomalous correlations 
play a crucial role on the stability of BEC and on the oc¬ 
currence of the superfluidity [15 ■ We find that the single 
bright soliton decays at nonzero temperatures, with the 
decay rate increasing with rising temperature owing to 
the host-host and impurity-host interactions. In addi¬ 
tion, we show that bright soliton trains can be produced 
automatically due to the modulational instability (MI) of 
the evolving classical phase in the impurity component 
of a harmonically trapped BEC-impurity mixture even 
with repulsive impurity-boson and boson-boson interac¬ 
tions. These trains generate without changing the trap 
geometry as has been suggested in [5 , neither using the 
Feshbach resonance as it has been observed in 3,11, or 
even without imprinting the initial wave function with a 
fluctuating phase as is shown in Q. The MI pattern as¬ 
sociated with the attractive interaction plays a key role 
in the formation of bright solitons in a pure BEC [15 [15 ■ 
By investigating the time evolution of soliton trains, we 
find that the number of bright solitons is increased with 
increasing the impurity-boson interactions. 

This paper is organized as follows. In section El we 
briefly review the main features of our theoretical ap¬ 
proach. Section uni is dedicated to analyze the behavior 
of a single soliton in the impurity component of a BEC- 
impurity mixture where we solve analytically and numer¬ 
ically the generalized self-focussing nonlinear Schrbdinger 
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equation. Section |TV] presents the generation strategy of 
bright soliton trains in a trapped BEC-impurity mixture. 
We show in particular how the solitary wave formation 
occurs in the impurity through the MI. Section |V] is de¬ 
voted to conclusion. 


II. FORMALISM 

We consider a mobile impurity of mass mj immersed 
in a BEC of atoms of mass ms at finite temperature. 
The impurity-boson interaction and boson-boson interac¬ 
tions have been approximated by the contact potentials 
gB5{r — v') and g/_B(5(r —r'), respectively. We neglect the 
mutual interactions of impurity atoms since we assume 
that their number and local density remains sufficiently 
small [l^ . fl^ and hence there is no impurity fluctuation. 
The TDHFB equations which govern the dynamics of 
the condensate, the thermal cloud, the anomalous den¬ 
sity and the impurity read 


ih^B = -b gB ((/3 - 2)nB + 2n + qn/)] (la) 

iUfh = 4 + gB {2Gm -I- 2n -b qn/)] m, (lb) 

ih^i = [h'"/+ gjBiuB+n)]<^i. (Ic) 

In the set (P), /i^ = —(?i^/2mB)A -b Vb and = 
— (?i^/2m/)A -b Vi are, respectively the single particle 
Hamiltonian for the condensate and the impurity, with 
Vb and Vj being, respectively the condensate and the im¬ 
purity trap potentials. <i>B(r) = (^/^^(r)) is the conden¬ 
sate wave function, ^^(r) = |$B(r)|^ is the condensed 
density, 4>/(r) = (^j(r)) is the impurity wave func¬ 
tion, n/(r) = |$ 7 (r)|^ is the density of impurity atoms, 
the noncondensed density h(r) and the anomalous den¬ 
sity m(r) are identified respectively with {r)'^{r)) — 
$g(r)$B(r) and {'ip{r)ip{r)) — $B(r)4>B(r), where 
and Ip are the boson destruction and creation field oper¬ 
ators, respectively. The total density in BEC is defined 
by n = n_B -bh. The dimensionless parameters /3 = U/gB 
with U = qB(l-bm/$|;) being the renormalized coupling 
constant G = PIA{[3 — 1) and q = giB/gB is the 

relative coupling strength. For /3 = 1, i.e., = 0, 

Eq. (fTal) reduces to the HFB-Popov equation which is safe 
from all ultraviolet and infrared divergences and thus 
provides a gapless spectrum. For 0 < /3 < 1, G is nega¬ 
tive and hence, rh has a negative sign. For /3 > 1, G is 
positive, and thus, rh becomes a positive quantity. 
Neglecting the mean-field interaction energy between 
bosons and impurity components i.e. giB = 0 and setting 
h = 771 = 0, one recovers the well known Gross-Pitaevskii 
equation describing a degenerate Bose gas at zero tem¬ 
perature and the Schrodinger equation describing a non¬ 
interacting impurity system. For further computational 
details, see Refs. l3ll - [^ . 

In our formalism the noncondensed and the anomalous 
densities are not independent. By deriving an explicit 


relation ship between them, it is possible to eliminate h 
via M 111: 


7= (2h-bl)^-4|r 


( 2 ) 


One can easily check by direct substitution that once 
Eq. ([5]) holds initially, it remains true during the dynami¬ 
cal evolution. The simplified set of equations are then the 
coupled equations (H)) with h is replaced by the expression 
(P|). In the uniform case, by working in the momentum 
space, Ik = coil? [ek/T) [^, where Sk is the excitation 
energy of BEC. At zero temperature 7 = 1 23 , and 
hence, Eq.Q reduces to h(h -b I) = Therefore, 

the expression of 7 clearly shows that rh is larger than 
h at low temperature, so the omission of the anomalous 
density in this situation is principally unjustified approx¬ 
imation and wrong from the mathematical point of view. 
Importantly, the expression of 7 not only renders the set 
© close but also enables us to reduce the number of 
equation making the numerical simulation easier. 

Moreover, what is important in the TDHFB approach 
for Bose systems is that there have been no assumptions 
on weak interactions. Therefore, the theory is valid even 
for strong interactions [3l|, . In addition, the TDHFB 

equations satisfy the total number of particles and 
the energy conservation law, and they provide a gapless 
spectrum 

The stationary TDHFB equations m can be eas¬ 
ily obtained within the transformations: ^B{x-,t) = 
^Bix) oqy^—igBt/h) and rh{x, t) = rh{x) exp{—igLmt/h), 
where g,B and /i* are, respectively, chemical potentials 
of the condensate and the anomalous density. If the con¬ 
densed and the anomalous densities change smoothly we 
can apply the Thomas-Fermi approximation, i.e. neglect 
the kinetic terms in Eqs. (fTal) and m- In this regime 
we immediately find that 


riB = 


(/3-2) 


m = 


2G 


9B 


-Vb „ 

- Zn — 777 ,/ , 

(3) 

9B J 

TA „ 1 

- Zn — jni . 

(4) 


The density profiles m and o have the form of an 
inverted parabola. For q = 0, an rh reduce to their 
usual expressions. It is remarkable that for ,5 = 1, the 
anomalous density vanishes. 


Assuming now that the impurity-boson coupling is suf¬ 
ficiently weak. Substituting Eqs. and (HI in Ea.dTB. 
we then obtain the generalized self-focussing nonlinear 
Schrodinger equation (NLSE) 


ih^i 



A + Vi 


XVb 


VI - gBlXni^ $/, 

(5) 


where A = q[I/(/3 — 2) -b 1/2G] and vi = 2ngBX — 

^[piB/W-2)-pirn/2G]. 
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The effective interaction potential A is positive only when 
the system is strongly correlated i.e for /3 > 2 or m > ub- 
Whereas it is negative for 1 < /? < 2. Importantly, A has 
a minimum at /3 = /3c = 6.83. 

Equation ([5]) has never been obtained before in the litera¬ 
ture. It allows us to study, in useful manner, the dynam¬ 
ical properties of the impurity and the fluctuation effects 
of the impurity environment at the same time. For /3 = 1, 
Eq.(l5]) reduces to the standard NLSE. 

Let us consider now an impurity immersed in an elon¬ 
gated (along the x-direction) BEC and confined in a 
highly anisotropic trap (such that the longitudinal and 
transverse trapping frequencies are iobxI^b^ <C 1). In 
such a case, the system can be considered as quasi-lD 
and, hence, the coupling constants of equations ([I|) effec¬ 
tively take their ID form, namely gB = 2hujB±aB and 
giB = ‘2,hu!B±aiB, where qb and qib are the scatter¬ 
ing lengths describing the low energy boson-boson and 
impurity-boson scattering processes. 

To well investigate the behavior of solitons in the im¬ 
purity component, it is convenient to reformulate Eq.® 
in terms of dimensionless quantities using the following 
parameters: x = x/lj where Ij = -^h/rrijuji is the im¬ 
purity oscillator length, a = mB/rni is the ratio mass, 
D_l = ojb±Iu}j^_, with w/x being the transverse impurity 
confinement frequency, r = tuji± and $/ = . 

Equation turns out to be given as 

-VI- gn^ $/, ( 6 ) 


where = 1 — AaD^, vi = vjjhuiix and g = 
2XV,±{aiB/li)- 

The energy functional corresponding to the NLSE m 
reads 


E = 






dx. 

(7) 


Remarkably, for strong pair correlations where /3 > 2, 
5 > 0, the impurity wave packet has an attractive 
interaction between atoms, so that the system supports 
automatically bright solitons. For 1 < /3 < 2, the 
impurity behaves as a dark soliton. 


III. SINGLE SOLITONS 

In this section, we restrict ourselves to study the prop¬ 
erties of single bright solitons. 

In the static case and for p = 0, Eq.® is integrable 
and admits the solution 

where A = 2/g is the extended localization length. The 
energy of the bright soliton is calculated through Eq. o 
as E = —(7^/24. 



FIG. 1. (Color online) Density profile of bright soliton in 
the impurity component. The parameters are set to: Nj = 
1 of ®®Rb impurity atom, of ^®Na bosonic atoms, 

a_B=3.4nm, ajs = 16.7nm, and the transverse trapping fre¬ 
quencies are ujb± = (jJix = 27r x 500 Hz. Blue dotted line: 
/3 = Pc, red dashed: line /3=3.5 and solid line:/3=2.5. 


To analyze the time evolution of the single bright soli¬ 
ton in detail, we fully solve Eq. (|6]) numerically employing 
the finite-difference splitting scheme. Figure. (P) de¬ 
picts that the soliton density decreases continuously un¬ 
til it reaches its minimal value for ,5 ~ /3c or equivalently 
TO = (/3c — l)ns ~ P-SPub- This latter condition which 
determines the soliton stability threshold, clearly indi¬ 
cates that the condensate is dominated by its anomalous 
fluctuations. As a result, a stable and robust impurity 
bright soliton locates at the local minimum of the effec¬ 
tive potential created by the impurity-host interaction. 
The decay arises from the fact that at nonzero tempera¬ 
tures the condensate coexists with both a noncondensed 
cloud and anomalous density composed of thermally ex¬ 
cited quasiparticles. Therefore, interactions between con¬ 
densed and noncondensed atoms on the one hand and in¬ 
teractions of the impurity soliton with atoms of the sur¬ 
rounding bath on the other hand lead to dissipation, and 
hence the soliton loses energy. Indeed, the effect of the 
interaction parameter /3 has considerable consequences, 
not only on the structure of the soliton, but also on the 
dynamics of the system, as we will see later on. 


IV. SOLITON TRAIN 

We investigate in this section the formation and the 
dynamics of bright soliton trains in the impurity compo¬ 
nent of harmonically trapped BEC-impurity mixture. 

Let us start first by illustrating the effect of the MI 
on the nonlinear evolution of the impurity. For a given 
stationary solution $/o the NLSE ([S]) with eigenvalue gi, 
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FIG. 2. The real and imaginary parts of the unstable eigen¬ 
value pertaining to the bright soliton versus the wave number 
of the corresponding eigenmode calculated by diagonalizing 
Eq. ([9)) (red circles). The blue curves correspond to the uni¬ 
form case as given by Eq. (nni). Parameters are the same as 
in Fig. [U 


the small-amplitude excitations of the system are defined 
through the random phase approximation (RPA) as $/ = 
<i)/o -I- 5$/ with (5$/ = -b 

We then obtain 


C M 
M C 


(Uq{x)\ 
\Vq{x)J 


= UJq 


( Uq{x) \ 
\-Vq{x)J ’ 


(9) 



X 



X 



X 


where C = —A/2 + p^a;^/2 — p. — 2gni with p = i? + fij 
being the chemical potential of impurity soliton, M = 
—and Uq(x),Vq(x) are the quasi-particle amplitudes. 
In the case of modulational instability analysis of homo¬ 
geneous quasi-ID case, the coupled system m is solved 
with being proportional to a traveling plane wave solu¬ 
tion of wave number k and complex dispersion ujq. Thus 
resulting in the following dispersion relation 


( 10 , 

For 5 < 0, the value of ujq is always positive, whereas 
for g > 0, the frequency ujq has an imaginary part 
and the spatially modulated perturbations grow expo¬ 
nentially with time, as is evident from the form of 
Therefore, the MI occurs for repulsive coupling constant 
g unlike to the case of a pure BEG. This can be attributed 
to the effects of the anomalous correlations which are 
strong enough to change the sign of the interactions. 
The fastest growth occurs for the wave number qmax 
that gives a maximum of ujq. The extremum condition 
dujq/dq = 0 gives qmax = and the maximum rate 

of growth is q^ — qmax/V^- In the context of the non¬ 
linear optics, MI occurs also for any sign of nonlinear¬ 
ity and GVD when the electromagnetic field is polarized, 
the light propagation in isotropic Kerr media is described 
by two incoherently coupled NLSE which referred to as 
cross-phase-modulation, and thus leads to MI for any sign 
of the interaction strength [s^ . 

In the trapped case, the dispersion relation can be 
obtained by taking the Fourier transform of the wave 


FIG. 3. Axial density profile of the BEG made for /3 = pc- 
Parameters as in Fig[T] 

number associated with each eigenmode calculated by 
diagonalizing the set The dispersion spectrum uiq 
is displayed in FiglJJ Plotting the analogous curves for 
the uniform case, where Eq. (nni provides analytical ex¬ 
pression, we observe an excellent agreement for the real 
part of the spectra while a small deviation appears in the 
imaginary part of the spectra at larger q. 

An important consequence of the MI is the formation 
of trains of bright solitons in the impurity component 
of a trapped BEG-impurity mixture. To investigate the 
generation of these trains, we first use the above ground- 
state wave function ([5]) as an initial condition for the 
time-evolution of Eq. ®. The numerical code imple¬ 
ments again in cylindrical symmetry a finite-difference 
splitting method. Assuming that the system is in equi¬ 
librium i.e. p = Pc (to avoid dissipation effects) and 
increasing the strength of the coupling between bosons 
and impurity ajB- 

Figure. [3] shows that at r = 5.4, the impurity cloud of 
®^Rb-^^Na mixture breaks into several peaks which re¬ 
main spatially localized during the time evolution and 
thus forming soliton trains. The number of peaks grows 
with the scattering length ajB- The formation of these 
trains can be explained as due to the strong anomalous 
correlations which grow themselves with increasing in¬ 
teractions [13 . Note that the pairing correlations lead 
also to split the bright soliton into two partially coherent 
solitonic structures of opposite momenta in a pure BEG 
[ 13 . Eurthermore, the bright soliton trains seen in Figl3] 
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generate by virtue of the MI of the time-dependent wave 
function of the impurity, driven by imaginary Bogoliubov 
excitations as we have mentioned above. 


V. CONCLUSION 


By using the TDHFB theory we have shown that bright 
solitons can be generated in the impurity component of 
BEC-impurity mixture at finite temperature as a result 
of strong anomalous correlations. We have found that 
these solitons decay due to the impurity-medium and 
medium-medium interactions. We have shown that after 
a sudden increase of the impurity-boson strength a train 
of bright solitons is produced in the impurity component 
and this can be interpreted in terms of the MI of the time- 
dependent impurity wave function. In our case the MI 


occurs even for repulsive interactions which is a salient 
feature in BEC-impurity mixtures. Most recently, the 
formation of bright solitons with repulsive nonlinearity 
has been observed from the surface spin-wave propagat¬ 
ing in yttrium iron garnet thin film strips 37|. Moreover, 
we have pointed out that the anomalous correlations, 
which become more and more stronger with increasing 
a/s, is also an important ingredient for the generation of 
bright soliton trains. The number of solitons in the train 
is proportional to qjb- 
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